Optical Mobius Singularities 
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Mobius strips with one, two, three, and four, half-twists are shown to be generic features of three- 
dimensional (nonparaxial) elliptically polarized light. The geometry and topology of these unusual 
singularities is described and the multitude of winding numbers that characterize their structures is 
enumerated; probabilities for the appearance of different configurations are presented. 
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PACS numbers: 42.25.-p, 42.25.Dd, 42.25.Ja 

The singularities of optical fields are major determi- 
nants of the field structure that are of intense current 
interest [1 — 13]. Here we extend previous studies of po- 
larization singularities in three dimensional (3D) optical 
fields [14— 22], finding new, highly unusual structures. 

3D light, i.e. light composed of waves with a broad 
spread of noncoplanar propagation directions, is almost 
invariably elliptically polarized. The generic singular- 
ities of 3D optical (and other electromagnetic) ellipse 
fields are lines of circular polarization, C lines, and lines 
of linear polarization, L lines [19 — 22]. At the point 
where a C line (L line) pierces a plane, S, a C point (L 
point) appears. The principal axes of the polarization 
ellipse are the major axis, the minor axis, and the el- 
lipse normal. At a C point the ellipse degenerates into 
a circle, the C circle, and the major and minor axes are 
undefined (singular). At an L point the ellipse collapses 
to a line, the major axis of the ellipse, and the minor axis 
of, and the normal to, the ellipse become undefined. 

All three principal axes of the ellipses surrounding the 
singularity, however, remain well defined. The projec- 
tions onto S of the major and minor axes (minor axis 
and normal) of those ellipses whose centers lie in this 
plane rotate about the C point (L point) with character- 
istic winding number ±1/2 (±1) [19 — 22]. Throughout 
this report, for C points (L points) we take S to be the 
plane of the C circle (to be oriented perpendicular to the 
major axis). We call this plane the principal plane, and 
label it Eo- 

Previous studies of 3D light have concentrated on pro- 
jections onto the principal plane Eo of the polarization 
ellipses whose centers lie in this plane; here we study the 
full 3D arrangement of these ellipses. Choosing for sim- 
plicity concentric circular paths in Eo that surround a C 
or L point, we find that for ~ 80% of all C points both 
the major and minor axes of the ellipses on each path 
generate a Mobius strip with a single half twist, whereas 
for the remaining 20% the Mobius strip has three half 
twists. Mobius strips on paths surrounding L points are 
generated both by the minor axes of, and the normals to, 
the surrounding ellipses, and have two full twists. Pre- 
viously [14], we showed that the polarization ellipses sur- 
rounding ordinary (nonsingular) points generate Mobius 
strips with one full twist, so that together with the results 
reported here, the polarization ellipses in 3D electromag- 
netic fields generically produce Mobius strips with one 



half-, two half-, three half-, and four half-twists. Repre- 
sentative Mobius strips are shown in Fig. [1] 




FIG. 1: Optical Mobius singularities in a random ellipse field. 
For clarity, the Mobius strips formed by the polarization el- 
lipses whose centers lie on a small circle in Eo are opened and 
straightened. Ellipse axes are shown by the thin straight 
lines that spiral around the thick straight line formed by the 
ellipse centers; only the upper half of each axis is shown. The 
small white (black) circle marks the ellipse at the arbitrary 
"beginning" ("end") of the opened Mobius strip; when the 
strip is closed these two ellipses are adjacent. For an odd 
(even) number of half twists the white and black circles are 
on opposite sides (the same side) of the line of ellipse cen- 
ters. (a)-(d) Mobius strips formed by the major axes (a)-(c) 
or ellipse normals (d) of the polarization ellipses surrounding: 

(a) a C point - this (wrinkled) strip has a single half twist; 

(b) an ordinary (nonsingular) point - this strip has one full 
twist; (c) a C point that is different from the one in (a) - this 
strip has three half twists; (d) an L point - this strip has two 
full twists. The Mobius strip in (a) forms a right-handed 
screw and has topological index —1/2, the strips in (b) - (d) 
form left-handed screws and have indices +1, +3/2, and +2, 
respectively. In random fields right/left-handed strips occur 
with equal probabilities. Although a circular path is chosen 
here for simplicity, any simple path that encloses only the C 
or L point yields (a possibly highly wrinkled) Mobius strip 
with the same twist index as does the circle 

Simple physical models of Mobius strips with the same 
topology as those in Fig. [T]can be constructed from short 
thin rods and a long strip of somewhat flexible material. 
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For a model of a C point (L point) the rods represent say 
the major axis (ellipse normal). The "upper" half of each 
rod is painted black, the rods are oriented parallel to one 
another and attached at their centers to the flexible strip 
such that they are perpendicular to the long axis of the 
strip. Introducing the appropriate number of half twists 
generates a model of the desired figure in Fig. [TJ gluing 
the ends of the strip together completes the construction 
of a model of the Mobius strip that surrounds a C or an 
L point. 

By continuity, for C points (L points) the angle that 
say the major axis (ellipse normal) of the surrounding 
ellipses makes with the principal plane goes to zero as 
one approaches the central singularity. In the physical 
model described above this can be accomplished simply 
by squashing the Mobius strip flat. Such a deformation 
changes the geometry of the strip but preserves its topol- 
ogy; the deformations in real optical ellipse fields are, of 
course, rather more subtle. 

In optical ellipse fields there are no packing problems 
due to overlap of the axes of neighboring ellipses in the 
concentric, nested Mobius strips surrounding C and L 
points. The reason is that the polarization ellipse de- 
scribes the time variation of the oscillating electric field 
at a point, so that the ellipse has dimension zero, whereas 
the Mobius strip has dimension one, and is a line. 

The Mobius strips in Fig. [T] were measured in a sim- 
ulated 3D random optical field (speckle pattern), where 
the optical field, E, in complex representation, was gener- 
ated as a sum of 100 plane waves with random phases and 
randomly oriented transverse linear polarizations. Many 
waves, however, are not required for C and L lines and 
their Mobius strips; all these features are found in three 
beam optical lattices. 

The C and L lines used in Fig. [I] were obtained from 
the above simulation by tracking zeros of the C and 
L point discriminants Dc = a\ — 4a2, and Dl = a2] 
these discriminants are obtained from the characteristic 
equation A 2 + a\\ + a 2 = of the real coherency ma- 
trix My = Re(E*Ej);i,j = x,y,z [16]. The orienta- 
tion of the major axis a and the minor axis (3 of the 
ellipses surrounding a C or L point were calculated us- 
ing formulas due to Berry [15,18]: a. = Rc(E*\/E • E), 
j3 = Im(E* VE • E); the normal to the ellipse, 7, was 
calculated using 7 = Im(E*xE) [19]. 

The three primary topological indices that characterize 
a C point in a 3D field are: (i) The well known [19] 

( o ( C} 

winding number I a = la = ±1/2 of the projection 
onto the principal plane So of the major or minor axes of 
the surrounding ellipses; (ii) The winding number I\ = 
±1 (not discussed previously) of the projection onto So of 
the normals to the surrounding ellipses; and (Hi) the first 
of the two major findings of this report: the twist index 

i~a — T fj — ±1/2, ±3/2 of the Mobius strips generated 
by the major and minor axes of the surrounding ellipses. 
For C points, axis 7 does not generate a Mobius strip. 
The above three indices generate 2 x 2 x 4 = 16 different 



C points, an eight-fold enlargement of previous C point 
characterizations [19]. All 16 C points appear in our 

simulations; T a is illustrated in Fig. Q] 1^ and 1^ 
are illustrated in Fig. [2] 




FIG. 2: C point indices I { a °> and if. I^ 1 (if 1 ) is gener- 
ated by the projections onto the principal plane Eo of axis a 
(axis 7) of the ellipses surrounding the C point. Convention- 
ally, the centers of the surrounding ellipses are taken to be 
located on a small circle surrounding the C point, although 
any simple closed path produces the same index. Projections 
of the axes of the surrounding ellipses onto Eo are shown by 
short straight lines, ellipse centers by small black dots, and 
the C point by the central black circle. As an aid in follow- 
ing the rotation of the axes one end is decorated by a small 
white circle. The here counterclockwise direction in which 
the surrounding circle (path) is traversed is shown by a curved 
arrow, the starting point of the path by a small straight ar- 
row, (a) for the C point in Fig. QJa). The net rotation 
of the ellipse axes is in the same counterclockwise direction as 

(C) 

the path, so the sign of is positive. The axes rotate by 
180° during one 360° traverse of the path, so the magnitude 
of the winding number is 1/2. Adding in the sign, yields 
lL C) = +1/2. (b) lL C) for the C point in Fig. Die). The 
ellipses rotate by 180° in the clockwise (retrograde) direction 
and lL C) = -1/2. (c) 4 C) for the C point in Fig. |]Ta). The 

ellipses rotate by 360° in the same clockwise direction as the 

(C) ' (C) 
path and = +1. I-, — +1 also for the C point in Fig. 

[2c). (d) if for another ellipse located on the same C line 

as the ellipse in Fig. []Jc). Here /i; C ' = — 1. 

The complex field components E x , E y , and E z . in the 
immediate vicinity of any point in a 3D field can be ex- 
panded to first order (the order required here) as E m = 
K ™ + J2n=x, v ,z ( p mn + iQmn)u n , where m = x,y,z, and 
u x = x, u y = y, u z = z. Without loss of generality, for 
C points we take K x to be real, set K y = ia \K X \, where 
a = ±1, and set K z = and u z = z = 0. This describes 
a C point at the origin with the normal to the C circle ori- 
ented along the z-axis, and the surrounding ellipses lying 
in the principal plane Sq of the C point, here the xy- 
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plane. Similarly, for L points we take K z to be real, set 
K x = K y = and u z = z = 0. This describes an L point 
at the origin with its major axis oriented along the z-axis, 
and again places the surrounding ellipses in the £y-plane, 
the principal plane So of the L point. In both cases, in 
the principal plane So, the generally nonzero parameters 
P nz and Q nz make no contribution because their mul- 
tiplier, z, is zero. The field component E z is not zero, 
however, because it involves terms in x and y that are not 
zero. We note that all structures found in our simula- 
tions were reproduced using this expansion with suitable 
choices for the expansion parameters K, P, Q. Analyt- 
ical expressions can be obtained for the various indices 
in terms of these expansion parameters by extending the 
methods of [14]; these methods require a discussion that 
is too lengthy to be included here. 

Close facsimiles of the measured random field strips in 
Fig. [T] can be generated from simpler expansions. The 
Mobius strip in Fig. [Ua), f° r example, is reasonably well 
approximated by Ex = —5 + y(i — 1), Ey = —5i + x, 
Ez = —ix + y; that in Fig. QJb) by Ex = —5, Ey = 
i, Ez — —y; the one in Fig. QJc) by Ex = 10 — y, 
Ey = lOi — x, Ez — y; and the strip in Fig. []Jd) by 
Ex = i (2y — x), Ey — x + iy, Ez = 5. 

Preliminary statistics of the various topological indices 
were obtained by choosing random values for the nonzero 
K, the various P and Q, and the sign of a. A random 
number generator that produced a normal distribution 
with unit variance was used. The data listed below were 
obtained from 1.6 x 10 6 independent samples. 

Fig. [3] presents a modified Zipf plot for the three C 
point indices l\ and • 
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FIG. 3: C point Zipf plot. Configuration fraction (probability) 
vs. configuration rank for configurations ranked by decreasing 
probabilities. 

There are four primary topological indices that char- 
acterize an L point. The first two indices are the well 
known [19] winding number 1^' = I^ L ' = ±1, and the 

index — ±1 (not discussed previously). These in- 



dices are analogous to the corresponding C point indices 
1^ = Ip and l!y G '. The other two indices are the 
second major finding of this report: the two independent 
twist indices, ri — ±2, and r-y = ±2. The Mobius 
strip in Fig. [ljd) was generated by axis 7 and has index 
Ty = +2, whereas the strip generated by axis (3 for the 
same L point has index — —2. For L points, axis a 
does not generate a Mobius strip. 

These four L point indices generate 2 4 = 16 different 
L points, all of which appear in our simulations. Fig. 0] 
presents a modified Zipf plot that summarizes the differ- 
ent L point configurations and their relative probabilities. 
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FIG. 4: L point Zipf plot. 

There are four geometrically distinct types of L points. 
We label these L mm , L mc , L cm , and L cc , where the first 
subscript refers to the figure generated by axis (3 of the 
surrounding ellipses, the second to axis 7; m (c) implies 
that the figure is a Mobius strip ("cone"). The dif- 
ferences between Mobius strips and cones are shown in 
Fig. [5] Although it is relatively easy to distinguish 
visually between Mobius strips and cones, it is difficult 
to do so with complete reliability using a computer pro- 
gram. From visual examination of some 100 L points we 
observed that L mm points appear mostly, but not exclu- 
sively, for negative 1^ and L cc points appear mostly, but 

not exclusively, for positive il ; L mc and L cm points ap- 
pear less frequently than do L mm and L cc points. The 
data in Fig. 0] do not distinguish between topological 
equivalent but geometrically distinct L points: this fig- 
ure lists all L points with a given index configuration 
without regard to whether these points are L mm , L mc , 
L C m 1 or L cc . 

The topological indices discussed here are all of first or- 
der in the sense that the field modulations that give rise 
to them grow linearly with distance from the singularity; 
these indices should be observable in microwave [3, 4, 20] 
and in interferometric, phase-sensitive, nanoprobe opti- 
cal [23 — 27] experiments that are able to probe the three- 
dimensional structure of the field. 
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FIG. 5: L point structures. The ellipses in Eo that surround 
an L point can generate one of two geometrically distinct 
structures: (a),(b) a double-full-twist Mobius strip, or (c),(d) 
a section of a cone. Ellipse half-axes of the surrounding el- 
lipses are shown by thick straight lines, the centers of these 
ellipses by light gray circles;the L point itself (not shown) lies 
at the center of this circle. The end points of the ellipse axes 
are shown by filled white (black) circles if these endpoints lie 
above (below) the plane of ellipse centers. In (a) the ellipse 
endpoints wind around the circle of ellipse centers to form a 
pair of interlocking rings. This Mobius strip is the one shown 
in Fig. 1(d). As shown here in (b), an observer walking along 
the circle of ellipse centers sees the endpoints spiral around her 
path (central straight line); the winding number of this spiral 
is r = +2. In (c) the endpoint circle and the circle of ellipse 
centers do not interlock, and this structure is not a Mobius 
strip. An observing walking along the circle of ellipse centers 
sees a section of an inverted cone whose " bottom" edge is her 
path and whose "top" edge, formed by the ellipse endpoints, 
is fluted and deeply scalloped; these endpoints form the spiral 
in (d), which has winding number r = +2. Inverted cones, 
such as the one shown here, and cones in which the ellipse 
axes flare outward, instead of inward, from the circle of el- 
lipse centers are equally probable, as are positive (+2) and 
negative (—2) values for r. 

C and L lines and their associated C and L points cor- 
respond to degeneracies of symmetric second order ten- 
sors (matrices) [17,18]. Accordingly, they, and their 
associated Mobius strips, can be expected to occur, and 
should therefore be sought, in other physical line and vec- 
tor fields described by such tensors; candidates include 
liquid crystals, strain fields, and flow fields. 

I thank Prof. David A. Kessler for helpful comments 
and suggestions. 
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